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ABSTRACT 


The purpose of this thesis was to develop a technique by which the 
frequency response of a servomechanism might be rapidly evaluated by 
graphical methods. The approach consisted of developing several frequency 
response parameters related to the location of the poles and zeros of the 
closed loop transfer function of a servomechanism. The systems considered 
were the pure second order (one pair of complex conjugate poles), second 
order with one zero (one pair of complex conjugate poles and one real 
zero), and a pure third order (one pair of complex conjugate poles and 
one real pole). For the latter two systems, the complex poles and the 
frequency were redefined as referred values relative to the real pole 


Or zero. 


Loci of constant frequency response parameters were mapped onto a 
complex plane; or loci of constant imaginary referred pole locations were 
mapped onto an auxiliary complex plane, depending upon the specific 
system considered. Using the points obtained from the appropriate 
relationship, the frequency response, M and N contours, could be rapidly 
and easily sketched. The M and N contours could then be used as a basis 


for rapid system evaluation in the frequency domain. 


The authors wish to thank Professor George J. Thaler for bringing 
to their attention the thesis by Major Walter B. Patton, USMC, and 
LT William B. Abbott, III, USN, without which this thesis topic would 


probably not have been considered /1/. 
The authors are deeply indebted to Assistant Professor Robert D. Strum, 
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agement, and advice given so selflessly. 
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1. INTRODUCTION 


Each year the importance and utilization of servomechanisms and 
automatic control systems increase. It naturally follows that more 
papers are submitted pertaining to the design, analysis, synthesis, 
and compensation of servomechanisms and automatic control systems. 

Many of these papers and existing texts discuss details and laborious 
calculations which are certainly a prerequisite to analyze or synthesize 
the operation of a servomechanism. However, in many cases, a rapid 
approximation of system operation would be beneficial, and possibly 

only an approximation would be required upon which an intelligent 
estimate could be made on the performance of the system under study. 
There are only a few papers known to the authors in which the use of 
approximations are presented from which an estimate of system performance 
can be made, e.g., Floyd's method /2/. One such paper by Patton and 
Abbott /1/ presents a method for a rapid approximation of servomechanism 


transient response. 


A corollary should exist to the Patton-Abbott development in which 
a rapid approximation of the frequency response of a servo could be 
determined. The use of frequency response analysis is still an important 
criterion in measuring servo performance, although the transient response 
approach has become very popular and widely used in recent years, almost 
to the point of neglecting frequency response. There are many ways of 
determining frequency response from transient response, but lengthy 


calculations are necessary. 


It is the purpose of this paper to present a graphical means for 
determining a rapid approximation of servomechanism frequency response. 
This approach eliminates the use of the Bode diagram and Nichols’ 
chart. The M and N contours will be obtained directly from the closed 


loop tranafer function of the servo. 
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This development will parallel the Patton-Abbott approach so that 


from the same given data (closed 
it is possible to obtain rapidly 
paper and the frequency response 
discuss in sufficient detail the 


approach so that no reference to 


loop transfer function, factored form) 
the transient response from their 
from this paper. This paper will 
development of the frequency response 


the Patton-Abbott paper will be required. 





2. GENERAL DEVELOPMENT 


Since the time domain or transient response is used so widely in 
the initial phases of servo study, all developments in this paper 
originated with the closed loop transfer function expressed in Laplace 
notation (s-plane) with all initial conditions equal to zero. The 
transformed equation was then converted into frequency response notation. 
For this to be a valid conversion, certain basic definitions are in order, 
plus some additional restrictions and assumptions used in this devel- 
opment. The frequency response of a system was defined in part by 
considering only steady state conditions. Thus, all transient conditions 
have vanished. To convert the transfer function from Laplace notation 
into the frequency response equation, it was only necessary to make a 
substitution of the variable. This was based on the assumption that 
the input was a pure sinusoid and that the output was also a pure 
sinusoid since the steady state condition exists. Furthermore, this 
forces the restriction that the system must be linear. The complex 
variable s = 0 + jw then becomes the imaginary frequency variable 
since the real part equals zero (0 = 0). Therefore s = jW; and when 
this variable substitution was made in the closed loop transfer 


function, the frequency response equation resulted. 


2.1 Basic Definitions and Systems Characteristics 
At this point, a tabulation of the various definitions, restric- 


tions, and assumptions would be in order: 

a. The system was linear. 

b. The system was stable. 

c. The system was characterized by unity feedback, i.e., 
M= 1 at w= 0. 

d. The closed loop transfer function was known in its factored 
form. 

Based on the above, some general frequency response characteristics 


could be stated. 





2.2 M-Contour 
a. At zero frequency, M is unity. 
b. At zero frequency, the slope of the M contour is zero. 
c. As the frequency approaches infinity, M approaches zero. 
d. As the frequency approaches infinity, the slope of the M 
contour approaches zero. 
e, At the maximum or minimum value of M, the slope of the M 


contour is zero. 


2.3 N-Contour 

a. At zero frequency, N is zero. 

b. As the frequency approaches infinity, N approaches fond 
degrees. 


(m is the number of zeros and n is the number of poles of the closed 


loop transfer function) 


2.4 Parameters 
To display graphically the frequency response, parameters were 
selected to define specific points on the M and N contours. The 


selected parameters were defined as; 

a3 MM The value of M where the slope of the M contour was zero. 
(Does not include the value of M at zero frequency when M *# 1.0). 

2) ee The value of omega (W * 2mf£) where the slope of the M 
contour was zero. (Does not include Ww = Q). 

c. N_: The value of N at ae 


p 
d. M The bandwidth point, or when M ® .707. 


Qe, Wo: The bandwidth frequency (when M # .707). 

ne New: The value of N at the bandwidth frequency. 

g- M: The notation when M #® 1.0, except at Ww # Q. 
Be W: The value of ™ when M # 1.0, except at W = Q. 
1. Ny: The value of N at W, - 


he 





ies D>: The undamped natural angular frequency. 


n 
lee Mun? The value of M at the undamped natural angular frequency aa 
iA Non: The value of N at the undamped natural angular frequency (WH ). 


The above selected parameters defined the general shape of the M 
and N contours. However, it was discovered that with only these points, 
the M and N contours were in some cases too general; the contours had 
insufficient shape to be of significant value. Better approximation of 
the M and N contours were provided by the selection of additional M, N, 
and W points to "fill-in" the shape of the contours. It was discovered 
that, mathematically, the choice of as the basis for these additional 
points reduced the M and N equations to relatively simple relationships 
and insured that the additional points were located in the critical 
regions of the contours. The benefits gained by the additional points 
more than compensate for the slight additional labor by providing a much 
more accurate approximation to the M and N contours. The additional 
parameters selected were: 

m. w: An arbitrary frequency proportional to 2 (w xk at 
where k was a selected constant as discussed later. 


n. M: The value of M at W 


O. Ny: The value of N at W 


A sketch of the general M and N contours is shown in Fig. 2-1. 


2.5 Classes of Servomechanisms Studied 

From Fig. 2-1, it can be seen that it was relatively easy to 
construct an approximate frequency response if the above defined para- 
meters were known. This paper contains the above parameters presented in 
graphical form. This paper also contains the technique utilized in the 
parameter determination such that the frequency response can be rapidly 
approximated for specific classes of systems. The three systems covered 
in this paper were: 

I. Pure Second Order - The closed loop transfer function consisted 


of two complex poles only. 





———— wee —_— aw 


- 180° 


We 4 


Wn 


Fig. 2-1 General M and N Contours Showing Parameters Chosen 





Il. Second Order with One Zero - The closed loop transfer function 


consisted of one pair of complex poles and one real zero. 


TII. Pure Third Order - The closed loop transfer function consisted 


of one pair of complex poles and one real pole. 


26 Presentation of Parameters 

The graphical presentations for the pure second order system were, 
for the most part, loci of constant values of the various parameters 
plotted on the s-plane. For the other two systems covered in this 
paper, it was necessary to express the coordinates of the complex poles 
as a function of the magnitude of the real pole or zero. In accomplishing 
this, the angular frequency, W, also became a function of the real pole 
or zero location. Thus the magnitude and phase shift for the frequency 
response was shown to be a function of the complex pole location relative 
to the real pole or zero. The graphical presentations for the second 
order plus one zero and for the pure third order systems were plotted in 
two ways. First, by loci of constant values of the various parameters 
plotted on a referred complex plane, that is, with ordinate and abscissa 
being the complex pole location relative to the real pole or zero. 
Second, by loci of constant referred imaginary part of the complex pole 
plotted against referred real part of the complex pole as abscissa versus 


the chosen parameter as ordinate. 


ae Digital Computer Utilization 


The solutions of the equations that will be presented later when 
a specific system is discussed were obtained from a Control Data Corp- 
oration high-speed digital computer Model 1604, located at the U. S. Naval 
Postgraduate School, Monterey, California (Appendix A). Without the 
availability of a high speed digital computer, it is doubtful that 
this paper could have been developed without the utilization of an 


analog computer plus excessive labor. 





3 THE PURE SECOND ORDER SYSTEM 


The pure second order system treated in this section was pre- 
viously defined as a system whose closed loop transfer function consisted 
of only two complex conjugate poles. This system has been extensively 
treated in most servo textbooks, usually in terms of the damping ratio 
and the undamped natural frequency. Here, however, it will be treated 
by the rectangular location of the closed loop poles. The closed loop 


transfer function of the pure second order system was stated as /i/: 


fs) _ OG; an 

© 7 d y (se * 2 . he by (3 : 1) 
The poles of this system were plotted on the s-plane in Fig. 3-1. 
Converting Eq. (3.1) to the frequency response form, it became 

ic (- = ae. 

iA a = a wi (352) 


" 
: t pears a) | L; 


Saleem and |W Calculations 


Equation (3.2) rearranged became 


—_e_Sse__arre—a—— : i ze 
c Z rie a © M — 


from which 


LAA 6 


Me SS 
YS f es oe Seas aE enka pope cag poem peemanat even 4 IE 
- { ; \e Fo : = ' 








(3.4) 


Note that if the system gain were constant, then "a" and "b'' would be 
fixed and the magnitude would be a function only of system frequency, 


l.e., M = £(wW) 








Fig. 3-1 Poles of the Pure Second Order System (s-plane) 





3.2 M andw Determinegtion 
_}2) ae See aia aaa 


ae was obtained by finding the value of W where the slope of the 
M versus W curve waS equal to zero. For matnematical convenience, this 
was done by first squaring Eq. (3.4), then taking the first derivative 


and setting the numerator equal to zero. 
4 ny] 4 Stale eee. aes 
-—(a+b) (2 +b -w )-2w) + 8a iv | =O (3.5) 


Simplifying, this becomes 
ee — 2 
a \b ae (3. 6a) 
and w#=0Q0 (see Section 2.4b) (3.6b) 


Substituting Eq. (3.6a) in Eq. (3.4) gave the value of M at ae 
i= — oe Ge) 


Note that from Eq. (3.6), Bh would not occur whenever "a'' was less than 
"pb", i.e., when the real part of the pole was greater than the imaginary 


Bart. (Figs. 3-2 and 3-3) 


Seo. WwW Determination 
bw 


The equation for the bandwidth frequency was obtained by substituting 
the bandwidth magnitude (M = .707) into Eq. (3.4) and solving for the 


frequency. This resulted in 





a, a a 
Co, Vi b-— ) +1414 ne (3.8) 
yy Vy - 7 


(Fig. 3-4) 


AOE 





3.4 M and W Determination 
(on n 


— 


The undamped natural frequency was defined as 





wr eya +b- (3.9) 


This equation expressed ae dieectly in terms of “al rand. Ub ‘eaihie. o- 5) 


The substitution of Ww ieeowade (3.4 el ded 


MM, a i. (3.10) 


Equation (3.10)expressed Ble in terms of the closed loop poles of the 


System. (Fig. 3-6) 


a5 W, Determination 


As defined previously, W, was the frequency, other than zero, at 


il 
which M was unity. Letting M = 1.0 in Eq. (3.4) and solving for w 


yielded 


1 


Simplifying this expression by recalling 2 Ae 


resulted in Cw, a7, Chef (oe) 


(Fig. 3-7) 


346 M and oe Determination 
As stated in Section 2, intermediate values of M and N, in addition 

to those determined above, would result in a smoother plot approximating 

the M and N contours. The values of k chosen for this class of servo 


were .5 and 1.5. These particular values were selected as being those 


-|]]- 





points which would be most beneficial in smoothing the plot of the 
M and N contours, and because the M equation was the same for k = .5 


as for k = 1.5. im oec tion 2, 


W, OC kW, 


k 
or ele <7), fa 
ae yk, Wy = 1.5 Coy 
i a 
and de = Ik, U4, =) { 5 Ud 
(mie 3-7) 
The value of Wy and W were substituted into Eq. (3.4) and yielded 
1 Z 
4 
M,= wie Ree 
| ! t ———= rhe) 
1. 5) ee 


-_ ae « 


- - CRS oa, 
Where M M. ie recalling wy ts 0. Loci of constant values 


or M,. were plotted in relation to "a" and 'b". (Fig. 3-8) 


3.7 N Calculations 


reon Eq. (3:3) 

[~ fe ra 4 Fae 2D “4 f 
Ke (j ww) = i se ee: ee = Whee 

i ae aera +22 } (DD 


from which, 


N Ae es af”: { 7) | Za Li) f 
= AYCLAN | = | (3.14) 


& +b wee 


ae 





Equation (3.14) was then the basic N equation for the pure second order 


system. 


ee 





3.8 N Determination 
__ i eee 


a was previously calculated by Eq. (3.6). The substitution of 
Eq. (3.6) into Eq. (3.14) resulted in the solution of ne related 


directly to "a'' and "b"’, the closed loop poles of the system. 


a. J (3.15) 


Loci of constant values of Me were then determined for all possible 
eo and “bl! yvalwes. (Fig. 3-9) 

It is pertinent to note that the equations presented were in most 
cases not solved as stated. Algebraic manipulations were performed to 
make the equation compatible to digital calculations by writing the 
equations in terms of "b" as a function of "a'' and the particular 


parameter. 


5 Determination 
SS en Sa 


N was determined in a manner Similar to Ni The value of > in 
terms of "a" and "b'"’ was stated in Eq. (3.10). This relationship of 


W was substituted in Eq. (3.14) and resulted in 





~ “~ 
N= arctan 2a 
— 6 bye | ERS Sel, Se is See ee ee oe 
Adj ae ae ee ae 4 
pe a Glatie | (3.16) 


which simplified to 








or 
—_— = 0 
N 90 


Thus, nee for a pure second order system was a constant and equal to -90 


oe 





Dero N Determination 
bw 


The equation for N__ in terms of "a" and 'b" was obtained by 


bw 
Siiestituting Bq. (3.8) into Eq. (3.14). This resulted in 





; 
eee 5 G 4 
N, =-aretan{2e' ee cg chee 


is Lab ald arab | = 


pond OS ya 


See. 3-10) 


Simi N, Determination 


Equation (3.12) determined the values of ™ when M = 1.0 for non-zero 


frequencies. When the expression of W™,, determined from Eq. (3.12), 


i tg 


was substituted into Eq. (3.14), 


(3213) 





Solutions to Eq. (3.18) were plotted as constant N., loci for all values 


1 
of "a" and "b"'. (Fig. 3-11) 


Be 2 Ny Determination 


As stated above, the k values selected were .5 and 1.5 where 


w, S¥k W49 « In a manner similar to the determination of N, above, 


§ 
owe aha \ 
=~— “ictani 


N ; uf _ 


—. Ps = an 3 
| 2a vk a | 


from which Ny values were determined explicitly in terms of "a" and "'b" 


for both values of k. (Figs. 3-12 and 3-13) 


3.13 Use of the Curves 
At this point, all the selected parameters necessary to sketch the 
M and N contours have been determined from the location of the closed 


loop poles. With the values of the closed loop poles known and the 


sae 





curves plotted as described above, a rapid graphical approximation of 
the frequency response of the pure second order system can be constructed. 
For convenience, the following recommended procedure is presented for 
the use of the pure second order curves: 

a. Enter Fig. 3-2 with the complex pole coordinates -a + jb, and 
obtain Ae Enter Fig. 3-3 with -a + jb, and obtain we (Note that oe 


is the ordinate in Fig. 3-3.) 


b. Similarly obtain a (M = .707) from Fig. 3-4, We and Mon 


mom Pigs. 3-5 and 3-6 respectively. 


c. OQObtain Wi (M = 1.0) from Fig. 3-7, entering with the oF 


determined above. 


d. Obtain M, FOr JDOEN ky = .5 and ko = 1.5 from Fig. 3-8, entering 
with "a" and "'b''. Determine nee (W, =V 5 “Gi,)s 
d f 1 Pas ih \ 
an “ky ® | @, =\1.5 sey 


from Fig. 3-7, entering with the value of sre obtained from step (a) above. 


The M contour for the specified closed loop poles may then be 
sketched, having easily obtained the M and corresponding W for six known 
points and realizing that M = 0 and has a slope of zero at W = Q, and 


also that the slope is zero at o and at W = ®, 


The N contour points are obtained in the same manner by using 
Figs. 3-9 through 3-12 to obtain the values of N corresponding to the 
same frequencies found for the M contour above. The one exception is 
Non? which is always equal to -90° for this system. Using these N, 

W values and knowing that N = 0° and has zero slope for W = 0, and also 
that N = -180° for w = © for a pure second order system, the N contour 


can now be accurately sketched. 


Note that this paper does not present curves which will give any 


values at frequencies higher than W since it is unlikely that there 


bw’ 
would ever be any interest in the region beyond the bandwidth frequency. 





3 14 Example 
Given: Tie ees 1 2 


Solution: 


SPECIFIC DATA 


TABLE OF PARAMETER VALUES 





GENERAL DATA (for all pure second order systems) 





Fig 3-14 (page 17 ) compares the frequency response obtained from 
this rapid approximation method to the actual M and N contours obtained 


from a digital computer program of the general M and N equations. 


Note: Although a scale is shown on these pure second order curves for 
convenience, the ioci are really independent of scale as they depend 
only on the ratio of "a" to ''b" Hence, any convenient scale for 9, 


j¥, and w may be substituted for the scales shown as long as all 


scales are in agreement. 
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4. THE SECOND ORDER SYSTEM WITH ONE ZERO 


The addition of a real zero to a pure second order system 


resulted in a system whose closed loop transfer function was 


“> b e . im 
fib (s td) 
F (sj2 BERS) (4.1) 


The singularities of this system were plotted on the s-plane in 


Fig. 4-1. The constant a + a , usually interpreted as System 
d 


gain, reflected the assumption of unity feedback as stated previously. 


By letting s = jW, the frequency response form of Eq. (4.1) 


resulted in 


- (ius) Se) ue! : -@- < Z a) 
al Pol Gea : ie” cA 4 . 
ed Gat ay eb (4.2) 





(dad b)43(wa+wb) 








e (4. 3) 
. (dar ade b~ dw ) +3 (2daw) | 
Dividing the numerator and denominator of Eq. (4.3) by a” 
Fs ie 2. , 2 
0- 91, c°f/a&a,~w dD 
(+ 3; +3 + ee) 
oe tat gril d* (4.4) 
(oe a o , ri | eo rn 
3 (2 4b _ wv) j (2 | .. w\ 





f 


From Eq. (4.4), it was seen that ‘a’, “b,and W may all be referred to the 


zero location, d. 
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je 


s-plane 


o~ 





Fig. 4-] Poles and Zero of the Second Order System with One Zero 
(s=plane) 
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From Eq. (4.4), 


[dl (4.5) 


b 
Mie Td] (4.6) 


a Ww 
ane a \d| (4.7) 


These quantities were illustrated in Fig. 4-2. Equation (4.4) then 


became 


feGod= 


(4.8) 





Thus the frequency response for the second order system with one zero 
became a function of the "referred" location of the complex roots and 
the "referred" frequency. This means then that the frequency response 
was the same for all systems having the same relative pole-zero con- 


figuration when the frequency was also referred to the zero location. 
(Fig. 4-3) 


4.1 M and Q Calculations 


Equation (4.8) was expressed in terms of its magnitude and phase 


angle by 


., _ Ln, a Me "1 + i (om, +m, a) 4 oh 
ca? Cm + Ma, od nF) +5 (2 yn 2) Me (4.9) 


ae 





s~plane 


1 | df 
idl 
~~ JM 
AG 3 
m-plane 
ee 
all 
Fig. 4-2 Referred Quantities of m, and m,> Second Order System 


with One Zero 
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system A: System B: 


a =@.0 a= .0l 
b = 2.0 b= .02 
d = 1.0 d= .01 
_ oe eS 
m, /a/ 1.0 m, /a/ lee, 
27D oe 
m5 /a/ 2.0 Mm, /a/ 2.0 





WW 
Salt 7) 
1 2 3 4 Ww, (d = 1.0) 
01 02 03 .04 w, (d = .01) 


Fig. 4-3 Frequency Response of Two Systems with the Same Relative 
Pole-Zero Configuration, Second Order System with One Zero 


age 





From this, 


M = V (m, ~ m2) +2 (m? m2) aie 
——»»» EEE a 
Von? +m, 0)” $4 


Equation (4.10) was the basic M equation for the second order system 


with one zero. 


4.2 M and Determination 
2) Se eee eee 


The equation for Me was obtained by first squaring Eq. (4.10), 
then differentiating with respect to {2 and setting the numerator of 


this equal to zero. This reduced to 


Ba ct V (m2 4m2) — 2 Cm? —-mz) + j 


Substituting this into Eq. (4.10) and simplifying, 


(ans) 





a 





Pp 


ee ee (4) 
me "+2, Mm, +m 2m, +2m2+1 +2 (nt, 1) 


Loci of constant values of m, were solved by digital computer and were 
plotted against m, versus Ms and i: (Figs. 4-4 and 4-5). 
4.3 Oe Determination 


The equation for the referred bandwith frequency was obtained by 
substituting M = .707 into Eq. (4.10) and solving for Ow After 


considerable algebra, this reduced to 
2 
_ 2 2 ee vs Q\2 2 2 
Uy 
1 2, \% a Z 
2 2. 
+ Gm toma) 1 |-4 moni f } (4.13) 


(Fig. 4-6) 
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4.4 M.~ and Qe Determination 


Qn 


The referred undamped natural frequency was defined for this 





system as 


4 
of) an 


ng = ee x (4.14) 


This expressed OQ as a function of m, andm,. (Fig. 4-7) 


Substituting Eq. (4.14) into Eq. (4.10) and simplifying yielded 


(4.15) 





(Fig. 4-8) 


a 5 2) Determination 


As defined previously, ©, was the referred frequency other than 


1 
zero at which M is unity. The equation for oa was found by letting 


M= 1.0 in Eq. (4.10) and solving for O 





= os ave 
AL, = | Gap etme) —-72 Oni — Mm; ) (4.16) 


(Fig. 4-9) 


4.6 M and 7 Determination 


The use of the curves mentioned up to this point for this class 
of system provided the four main points required to Sketch roughly the 
M contour. In order to "fill-in" the blank portions of the contour on 
both sides of ae it was decided to provide curves for three additional 
points on the M contour in the same manner as was done for the pure 
second order system. The use of these additional points was not 
necessary to sketch the frequency response; however, their use provided 


a more accurate approximation of the M contour. 


2362 


ae 





As in the pure second order system, the following definition was made 


= Vk $1.5 (4.17) 


wnere k is a constant. 


SUbDSEItULInNg into Eq. (4.11) for ae 


i= b 
% 


lot gaia 2 ae ee ae te 
{), = | -% Tk | (Mt mg | “LIM, mm, Teo (4.18) 


Values of .5, 1.5 and 2.5 were selected for k. M solutions were 


obtained for each k value by solving Eq. (4.18) for a and then sub- 
Stituting this value into the basic M equation with the same values of 


m, and m,- This was accomplished handily by one digital program. 


(Figs. 4-10, 4-11, 4-12,and 4-13) 


4.7 N Calculations 
Referring to the basic closed loop transfer function in frequency 
form, Eq. (4.9), the general phase angle equation for the second order 


system with one zero, was written as 


—_ 


N= arvetan see (4.19) 


+ 2 ow L My age 


The previously obtained parameter frequency equations Oo a a 


0,» and the three ys) were solved by digital computer programs. The 
resulting parameter frequencies were then substituted into the general 


N equation (4.19) along with the corresponding values of mM) and m.,- The 


following figures resulted, all plotted with constant mM. loci against m, 


versus N: 


Fig. 4-14 Ny as function of mM, and m., 


Fig. G15 Naw ais funetionvor my and m., 


Fig. 4-16 No, as function of Mm and m., 


Fig. 4-17 N. as function of m, and m 


1 1 2 

Fig. 4-18 N as {uncllon. Of mm. andem tone k se — ee 
ky} 1 2 1 

Fig. 4-19 N, as function of m, and m, for k, = 1.5 
k9 1 Z 2 

Fig. 4-20 N as functien-ef mand mm. Pore 2 


k 3 1 2 3 


wey) = 





4.8 Example 


Givemm™ Complex poles; -a + jb = -2.0 + j 2.0 


Real zero; -<-d =~2.0 


Solution: 


SPECIFIC DATA 


TABLE OF PARAMETER VALUES 


P 
bw os 
| 142 
2.0 
ky 2.78 
k. 
2 || eee 


“3/274 





---~ 
Obtain W by w = /d/Q 
x x x 


GENERAL DATA (for all Second Order plus One Zero Systems) 





Fig 4-21 (page 39 ) compares the frequency response obtained from 


this rapid approximation method to the actual M and N contours 


obtained from a digitai computer solution of the general M and N 


equations. 
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») THE PURE THIRD ORDER SYSTEM 


The addition of a real root to the pure second order system 
resulted in a new system which will be referred to as the pure third 
order system. The closed loop transfer function of the pure third 


order system had one pair of complex conjugate poles and one real pole. 


The closed loop transfer function for the unity feedback third 


order system was written as 


et bie 


[> (shox —S ee a; : 
cee fsa cetadibn seal) oa 


Fig. 5-1 defined these roots on the s-plane. 


Equation (5.1) was simplified as follows: 








P j ae 
[. (s)= === oe 
¢: : ‘ A a . ( : 
yo pts f 4 4 ” ot b f 
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~~ % ie | See if & Peat 4 San 2 2 
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Letting s = jw in Eq. (5.3) resulted in the frequency response form 


of the transfer function 
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Numerator and denominator of Eq. (5.5) divided by ne 


bese 
i (=) ca aa = | X\ ean (5.6) 
“s < jl a j fae i ne : 4 ‘ e a : P P a A at e 
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pre-e tS OS |) ONG BV Cee ce 


we ou 
Observing that a power of c exists in all terms corresponding to 


the same power of "a", "b'' and w; the following definitions were 


made: 
ao 
n=— Dick 
1 le] Or) 
A bh 
ie (5.8) 
wc! 
A wo 
= a) (5.9) 
These quantities were illustrated in Fig. 5-2. Equation (5.6) then 
became 
7k) ie ee A iN 
LGw= aa * ms ean ; —_ oan Pt _ Wire 
. ee ae ea Me ee - 
[ re, +12,- 9. (210,41) | +5. Gi74 m, +72, J 2 | (5.10) 


ny and n, were the complex pole locations referred to the real pole 
location and &2 was the referred frequency. The frequency response 
for the pure third order system was then a function of the referred 
complex pole and the referred frequency. Thus all systems having the 
Same relative pole configuration will have identical response if the 


frequency were also referred to the real pole location. 


5.1 M and © Calculations 
The magnitude portion of Eq. (5.10), after algebraic manipulation, 


reduced to 
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Fig. 5-2 Referred Quantities ny and N»» Pure Third Order System 
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Equation (5.11) was the basic M equation for the pure third order 


System. 


5.2 Discussion of the M Contours 

It was found that the M contour for a pure third order system 
will be one of three basic "shapes" as shown in Fig. 5-3a, b and c, 
depending on the values of ny and n,- Figure 5-3a was the situation 
when there was no peak. No value of ee was obtained and hence no values 
for oF or the (os. Figure 5-3b was similar to the second order M 
contour shape and had a single peak. For this general "shape" points 


e ie , ; 
were obtained for Le 0? ie Ca, Heo? and Ae 


The third general ''shape" is Fig. 5-3c, which resulted from the 
. equation (5.12) which yielded two non-zero frequencies at which 
the M contour slope was zero, one being a "peak" frequency and the 
other being a ''minimum'' frequency. This latter frequency was designated 
Qe to avoid being confused with the oa It was also discovered that 
with this shape of M contour, it was possible to have 2 less than 
unity. Another possibility was to have Mh less than Mow and at the same 


time MS greater than one. This resulted in two values of & two values 


Le 
C +f 12 e 

of ky’ Me? ee and nie (one relative to M, and the other relative to 

M > and three values of vee (one for each crossing of the M = .707 line), 


in addition to the usual points at ne Ca and QQ. 


Figure 5-4 might be thought of as the complex n-plane and shows the 
regions in which certain ny and nN» combinations produce one of the 
three above M contours. This figure indicates the general shape of the 
M contours that result from the know ny and ny values. 


5.3 M and Determination 
2 Se ee ee 


The equation for was obtained by squaring Eq. (5.11), then 
differentiating with respect to {2 and setting the numerator of this 


equal to zero. This reduced to 





(5.12) 
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Fig. 5-3 The Three Basic M Curve Shapes for Pure Third Order System 
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The plus and mipus values of the inner square root quantity were both 
valid and resulted in two possible solutions for The consequences 


of two solutions were discussed in detail in paragraph 5.2 above. 


Loci of constant values of n, were solved by a digital computer 


program and were plotted against n, versus se and oa in Figs. 5-5a 


1 
and 5-5b In the same computer program, values of Mt were obtained 


for corresponding values of ny» ee and n, by substituting these 


known quantities into the basic M equation, Eq. (5.11). Figures 5-6a, 


5-6b, and 5-6c contain a plotted as a function of ny and n,- No 


specitic Me equation was derived due to the obvious complexity 


mivaived in substituting Eq. (5.12) into Eq. (5.11) for %&%. 


Bea 8 «C, Determination 
“bw 


The equation for the referred pandwidth frequency was obtained 


by substituting M = .707 into Eq. (5.11) and solving for n, as a function 


of Oo and ny» This involved extensive algebra which reduced to 


a Pad 
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op ff Mbt My A412) 4 Fig ty LL) = 
Le a Lind 
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wy i 


The + signs were both valid here also and their effect will be discussed 
Hater. (Fig. 5-7) 


bm  Ms~-and O° Determination 


On 
The undamped natural frequency for the pure third order system was 
defined as 


—— a ee ee Set) 
. ‘ 4 


dy J 
AX 
D at i ow 


Ny, 2 Ve a 


OQ was expressed as a function of n, and n,. (Fig. Sesh) 
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Substituting this into Eq. (5.11) and simplifying yielded 


en 2 a 
a f sou 
ne 


Df Fe RONG ts eg : 
V 42 (n, +n Mn, ser Mies 1) (Geis) 
Mon was thus expressed as a function of ny and No: (Figs. 5-9a and 5-9b) 


5.0 SF Determination 


The referred frequency, other than zero, at which M ® 1.0 was 
defined as 0. 


Letting M = 1.0 in Eq. (5.11) and solving for 0 resulted in 





a ok a 
N,~7, -4N nN, +.25 (5.16) 


(Figs. 5-10 a, b, c, and d) 


Do. 7 M. and om Determination 


For the same reasons as discussed in Section 2.4, it was found 
desirable to obtain "fill-in'' points on both sides of vo This 
was done by defining 
; QT = kOQp (5.17) 


which eliminated the square root of k as used for the previous two 


classes of systems. 


The values of k chosen were .8, 1.4, and 1.8; therefore 


= .8 Mp (5.18) 
ey 

#1.4 Mp (5.19) 
Mey 
a, = 1.8 Gp (5.20) 


G4 





eo egret La —— 


Fig. 5-ll was the almost trivial plot of these three equations, 
included for optional use in determining the cy values once = has 
been determined. 

Loci for the three M's were calculated by a digital program which 
Meta zed Eqs. (5.11), (5.12), (5.18), (5.19), and (5.20); thus it was 
not necessary to attempt to obtain any specific mM equations. 


(Figs. 5-12 a,b, 5-13 a,b, 5-14) 


5.8 Determination of N Curves 
Referring to the basic closed loop transfer function in frequency 
form, Eq. (5.10), the general phase angle equation for the pure third 


order system was written as 


N=-arctan{ =Een Cres n') se 
~" (21 }-j rt. ek 


The previously obtained parameter frequency equations oe OA 
a and a) were solved by digital computer programs. The resulting 
parameter frequencies were then substituted into the general N equation, 
Eq. (5.21), along with the corresponding values of ny and N,- The 


following figures resulted, all plotted with constant n, loci against 


2 


n, versus N: 


1 


Pag. o-15 a Me as a function of n, and ny 
Fig. 5-15 b N as a function of n, andn 
m 1 2 
Fig, 3-16 a Now as a function of ny and n, 
Big. o=16. b Naw (extended) as a function of ny and n, 
Fig. S-l7 a Ny as a function of ny and n, 
Bic. 5-17 b Ny (expanded) as a function of ny and ny 
hae. O- 1/7 ¢ Ny (2nd valuejas a function of ny and n, 


Fig. 3-18 a N, for maximum OE 
] 
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Ricseo-16 b N, for minimum (2 ) 
1 
Fig. 5-19 a She, for maximum oo 


Fig. 5-19 b Ny for minimum (QQ ) 
9 m 


Fig. 5-20 N. for maximum (£2 ) 
3 P 


An equation for Non was obtained by substituting Eq. (5.14) into 


the basic N equation, Eq. (5.21), 


f 
Sagi | eee We 
N, = arctan 4] 


LA 


| =a 
= —arctarn lel (5.22) 


i my | o> 2. 
re 13-4 “ te ayes Ped 

Vit, $7, 
ow 


oe1 of Non were plotted as a function of ny and Ny (Fig. 5-21) 


B66. 





5.9 Example 








Given Cemex POlCS: “Saat jo = 6 Tow 
Real pole: -c =-2.0 
solution: 
a 
SPECIFIC DATA Des ea 4 
a 20 
O77 ie een 


TABLE OF PARAMETER VALUES 





P £88 
mn CGS 
bw 2 34 
bw | NA. 
bw N.A. 
n 2.04 
: 14 
1 Veo 
ky es Z.76 
“ll Aba 
te 
3 a. 7 
* Relative to {. ** Relative to Q wae Ww = f/ef ©. 
m x x 
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GENERAL DATA (for 211 pure third order) 





Fig. 5-22 (page 69 ) compares the frequency response obtained from 
this rapid approximation method to the actual M and N contours 
obtained from a digital computer program of the general M and N 


equations. 
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6. CONCLUSIONS AND RECOMMENDATIONS 


The following conclusions were made based upon the investigation 
of the three classes of servomechanisms studied in this thesis: 

a. All pure second order servomechanisms with the same relative 
complex pole location will have the same frequency response. 

b. For the pure second order servomechanism, the parameter loci 
plotted are made applicable to any location of complex poles in the 
left half of the s-plane by changing the scales on the complex plane 
and on the parameter frequency loci. 

c. The frequency response for servomechanisms with three 
Singularities can be expressed as a function of the referred real and 
imaginary parts of the complex poles and the referred frequency where: 

1. The referred real part of the complex poles is the 
actual real part of the complex poles divided by the 
magnitude of the real pole or zero. 

2. The referred imaginary part of the complex poles is the 
positive imaginary part of the actual complex pole 
divided by the magnitude of the real pole or Zero. 

3. The referred frequency is the actual frequency divided 
by the magnitude of the real pole or zero. 

d. All servomechanisms with three singularities which have the 
same relative pole-zero locations will have the same frequency response 
relative to the referred frequency. 

e. Curves for determining values of several optional "fill-in" 
parameters were included in addition to the normally accepted parameters 


Ww | z 
of Mw? 5 and WD, ag 


f. An accurate graphical approximation of the frequency response 
of a servomechanism can be rapidly obtained uSing the curves produced 


in this thesis. 


It is recommended that the following areas be considered for further 
investigation using the technique presented in this paper; 
a. Servomechanisms with more than three singularities, and 


b. Systems other than servomechanisms. 
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8. Appendix A 


The Control Data Corporation Model 1604 digital computer 
important characteristics are summarized as follows: 
1. All transistorized 


2. 48 bit word length 





3. Stored program 


gS 


Single-address logic, two instructions per 48-bit word 


5. Indirect addressing feature 





| 
| 6 32,768 48-bit words of magnetic core storage: 


(a) Stored in two independent 16,384 word banks, 
alternate phased. 


(b) 4.8 micro-seconds effective cycle time 
(representative program). 


| (c) 6.4 microseconds total cycle time. 
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